Dynamics of thermal Casimir— Polder forces on polar molecules 
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We study the influence of thermal Casimir-Polder forces on the near-surface trapping of cold 
polar molecules, with emphasis on LiH and YbF near an Au surface at room temperature. We 
show that even for a molecule initially prepared in its electronic and rovibrational ground state, 
the Casimir-Polder force oscillates with the molecule-wall separation. The non-resonant force and 
the evanescent part of the resonant force almost exactly cancel at high temperature which results 
in a saturation of the (attractive) force in this limit. This implies that the Casimir-Polder force on 
a fully thermalised molecule can differ dramatically from that obtained using a naive perturbative 
expansion of the Lifshitz formula based on the molecular ground-state polarisability. A dynamical 
calculation reveals how the spatial oscillations die out on a typical time scale of several seconds as 
thermalisation of the molecule with its environment sets in. 

PACS numbers: 34.35. +a, 12.20.-m, 42.50.Ct, 42.50.Nn 
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I. INTRODUCTION 



Cold ensembles of polar molecules such as YbF have 
recently received particular attention due to their poten- 
tial use as ultra-sensitive probes of the permanent electric 
dipole moment of the electron [l[ , measurements of which 
allow for investigating the possible existence of physics 
beyond the standard model Q. The need for longer in- 
terrogation times has led to the development of Stark 
deceleration techniques for these heavy molecules [3, 0] , 
with a view to ultimately be able to trap molecules near 
microstructured surfaces (chips). Recently, trapping of 
light molecules such as metastable CO in travelling po- 
tential wells near a chip surface was achieved [f| . Another 
light diatomic molecule that has received considerable at- 
tention due to its large dipole moment is LiH; and the 
production of supersonic beams of cold LiH has been re- 
ported [6j. 

When attempting to trap polar molecules in close prox- 
imity to a surface, attractive Casimir-Polder (CP) forces 
0| ■ - effective electromagnetic forces between a neu- 
tral and polarisable particle and a macroscopic object 
- need to be taken into account as an important lim- 
iting factor. Thermal CP forces on atoms at thermal 
equilibrium with both the electromagnetic field and the 
present macroscopic bodies have been intensively stud- 
ied in the past on the basis of Lifshitz theory [8l-tllj , lin- 
ear resp onse theory [l2|, [HJ or normal-mode techniques 
[l4L [HI]. At room temperature, the energies associated 
with atomic transitions are much larger than the ther- 
mal energy, Hloa S> k^T, resulting in very low thermal 
photon numbers. A "high-temperature limit" is only ac- 
cessible in a geometric sense when the atom-surface sep- 
aration za is much larger than the thermal wavelength, 
za 3> he/ (2Tik-QT); in this case the thermal CP force on 



the atom can be approximated by [12l . \lc 
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for a two-level atom (transition frequency ujai dipole ma- 
trix element (1a) interacting with a perfectly conducting 
plate (unit normal e z ). 

The situation is different for molecules: whereas tran- 
sition energies of atoms are typically much larger than 
attainable thermal energies, the energies associated with 
rotational and vibrational transitions of molecules, heavy 
molecules in particular, are often small compared to the 
thermal energy even at room temperature. A genuine 
high-temperature limit Hloa *C T is hence realised with 
an associated large number of thermal photons being 
present. An additional consequence of the long transi- 
tion wavelengths is the fact that CP forces on molecules 
are expected to have a long range with the non-retarded 
regime za <C c/loa extending quite far out from the sur- 
face. A naive application of the above formula ((T|) for 
atoms beyond its scope to the high-temperature limit 
hh)A <C k-gT would suggest that the force can get arbi- 
trarily strong for molecules of smaller and smaller transi- 
tion energies, which already indicates that CP forces on 
molecules must be treated with care. 

Supersonic beam expansions typically produce cold 
molecules that are to a large fraction in their rovibra- 
tional ground states. For example, in the experiment re- 
ported in Ref. [(| , 90% of the observed cold LiH molecules 
were in their electronic and rovibrational ground state 
Y 1 S + . The cold molecule and the room-temperature 
surface are thus strongly out of equilibrium with respect 
to each other, so a study of the CP interaction necessi- 
tates that account be taken of the full non-equilibrium 
dynamics of the rotational and vibrational degrees of 
freedom of the cold molecule coupled to its thermal envi- 
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ronment. In contrast, in the context of non-equilibrium 
forces on thermalised atoms in an environment of non- 
uniform temperature, as recently proposed [l(| and mea- 
sured [13], a study of the full internal atomic dynamics 
was not necessary. 

In this article, we study the non-equilibrium thermal 
CP force on a polar molecule which is initially in its elec- 
tronic and rovibrational ground state in the vicinity of a 
metal surface. A recently developed dynamical theory of 
forces on single atoms or molecules in arbitrary internal 
states and arbitrary uniform temperature environments 
[l8| provides the necessary framework (note that a similar 
theory has been developed for two-atom van der Waals 
forces 19]). In particular, we will show that in contrast 
to the above intuitive expectation obtained from com- 
parison with the atom case, the attractive CP force on a 
molecule saturates in the high-temperature limit. 



II. CASIMIR-POLDER FORCE FOR GIVEN 
MOLECULAR STATES 

We consider a polar molecule (energy eigenstates 
eigenenergies huj n , transition frequencies u mn =ui m — ui n , 
dipole matrix elements d mn ) which is prepared in an 
incoherent superposition of its energy eigenstates with 
probabilities p n . As shown in Ref. the thermal CP 
force on such a molecule is given by 



n 

with perturbative force components 

oo 
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Here, G^ 1 ' is the scattering part of the classical Green ten- 
sor for the electromagnetic field in the given environment 
and £n = ^irk^TN/h denotes the Matsubara frequencies. 
The CP force ((3| contains both non-resonant contribu- 
tions (first term) and resonant ones (second term), where 
the former would also follow from applying Lifshitz the- 
ory in conjunction with the ground-state polarisability 
(we refer to it as Lifshitz-like force in the following) and 
the latter are due to the absorption and emission of ther- 
mal photons with photon number 



n(uj) 
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Given a probability distribution p n , Eq. ([3]) allows us 
to compute the thermal CP force. In particular, if the 
molecule is in an isotropic state such as the ground state 
or a thermal state (see below), the forces simplifies to 

m 
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A. Molecule near a plane surface 

For example, let us consider a molecule at a distance 
z A from the planar surface of a (non-magnetic) substrate. 
The respective scattering Green tensor is given by [2(| 



P 2 <? \ , , og 2 c 2 

2 J K^x^-x \ €y€y) ~\~ 2 TpG z G z 



where 
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e{u)P ~ Pi 
e(w)/3 + Pi 
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with Im/3, Im/?i > are the Fresnel reflection coeffi- 
cients for s- and p-polarised waves, /3 = y/uP-jc? — q 2 
and Pi = \J s(lS)u3 2 I c 1 — q 2 are the ^-components of the 
wave vectors in free space and inside the substrate, and 
e(u) is the (relative) permittivity of the substrate. Sub- 
stitution of G*- 1 ' (r, r, ui) into Eq. ([6| above leads to an 
explicit form for the CP force. 

The results simplify in the non-retarded and retarded 
limits of small and large atom-surface separations. In 
the non-retarded limit maxj{| y/ e(cJi) \u}i}z A / c <C 1 {ujf. 
relevant molecular and medium frequencies) , the approx- 
imation p ~ pi ~ iq leads to 

F n (r A ) = T > {l-^6 N0 )a n (i£ N ) jrHr—. e z 

|e(w„ fe )| 2 - 1 
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Note that while applying well to dielectrics, the non- 
retarded limit often provides a very poor approximation 
for metals, because the large factor |y / e r (wi)| may restrict 
its range of applicability to extremely small distances. 

In the retarded limit u> m i n ZA/c 3> 1 (wmin- minimum 
of the relevant molecular and medium frequencies), the 
resonant part of the force is well approximated by letting 
q ~ 0, while the approximations a n (i£jv) ~ a n (0) and 
e(i^Ar) ~ e(0) hold for those £n giving the main contri- 
bution to the non- resonant part. The q- integral for the 
N = 1 term can then be carried out immediately, while 
those for the remaining part of the sum can be rewritten 
in a more convenient form by introducing the integration 
variable v — /3c/£jv. Performing the sum according to 
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where x = 2irk B TzA/{lic). In particular, for a conduc- 
tor whose plasma frequency u>p is large compared to u> nk 
[cf. Eq. (|14p below] one has |e| 3> 1 and the retarded CP 
force is well approximated by 
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Note that the retarded limit as given above holds for 
all distances which are sufficiently large with respect to 
the atomic and medium wavelengths, irrespective of the 
temperature. When in addition the distance is very large 



with respect to the thermal wavelength (such that x ^> 
1), the contribution from the second terms in the above 
Eqs. (fT^l) and (ITSl vanishes and the non-resonant force 
approaches its well-known (geometric) high-temperature 
limit, cf. Eq. (JTJ) . In the opposite case of a distance which 
is much smaller than the thermal wavelength (x <C 1), 
the first terms vanish and the non-resonant force reduces 
to its (retarded) zero-temperature form, cf. Ref. [2l|. Our 
results, in particular those for the non-retarded limit, 
agree with the ones previously obtained in Ref. . Note 
that resonant force components and their oscillatory be- 
haviour in the retarded regime were first discussed for 
excited atoms at zero temperature (cf., e.g., [H, HI]). 

The limits reveal that the CP force follows a l/z\ 
power law for non-retarded distances. In the retarded 
regime, the non-resonant force components again follow 
an inverse power law whereas the resonant force compo- 
nents give rise to spatially oscillating forces whose am- 
plitude is proportional to 1/za- If present, the resonant 
force components are dominating over the non-resonant 
ones, in general. The magnitude of the contributions 
from various molecular transitions to the force ^ is de- 
termined by their dipole matrix elements and frequen- 
cies, where Eqs. (fl"U|) and JT2|) together with Eq. ([7]) im- 
ply that the strength of the non-resonant force is roughly 
proportional to l/uj kn , while that of the resonant force 
is governed by n{u nk ) + 1 or n{uj kn ) in the non-retarded 
limit and by w 3 fc [n(w„fc)-|-l] or ujf. n n(iu kn ) in the retarded 
regime. Equations pH|) and (fT!?j) furthermore show that 
the force becomes larger for larger permittivity of the 
surface material and saturates in the high-conductivity 
limit. 

The general results and discussion given above can be 
easily applied to various polar molecules interacting with 
different surface materials. The qualitative behaviour of 
the forces will be similar for all molecules and materi- 
als, i.e., a power-law dependence for non-retarded dis- 
tances will give way to an oscillating force in the re- 
tarded regime. The exact magnitude of the force as 
well as the length scale of the oscillations will depend 
on the dipole moments and frequencies associated with 
the specific molecular transitions involved and the elec- 
tric response of the surface in the way indicated above. 
Tabulated data for a variet y o f molecules and metal sur- 
faces can be found in Ref. [2J] . In the following, we will 
consider two representative examples. 



B. Examples: LiH and YbF near an Au surface 

We first consider a LiH molecule in its electronic, vi- 
brational, and rotational ground state (p n = 5 n o) near 
an Au surface at room temperature T = 300 K. With 
the help of the Green tensor ([5]) , we are able to compute 
the force components according to Eq. ([6| which are dis- 
played in Fig.[T] For this molecule, the contribution from 
rotational transitions with u> kn = 2.79 x 10 12 rad/s and 
d = 1.96 x 10~ 29 Cm [H {Y, k d ok d m = d 2 \ where I is 
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FIG. 1: Thermal CP force on a ground-state LiH molecule 
near an Au surface. For details, see text. 
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FIG. 2: Thermal CP force on a ground-state vs fully ther- 
malised YbF molecule near an Au surface. For details, see 
text. 



the unit tensor) is strongly dominant over those of vibra- 
tional and electronic transitions with their considerably 
higher transition frequencies. Molecules with a similar 
behaviour include NH, OH, OD, NaCs and KCs. For the 
relative permittivity of the Au surface we have used a 
Drude model 



e(w) = 1 



uj{uj + vy) 



(14) 



with w P = 1.37 x 10 16 rad/s and 7 = 5.32 x 10 13 rad/s 
[25] ]. In view of the current debate regarding the ther- 
mal Casimir force (cf. [26|, H3] and references therein), 
we have also calculated the force using the alternative 
plasma model and found that the difference between the 
two models is of no importance in our case. 

Figure Q] shows the contributions from the non-reso- 
nant force component (thin solid line) which is seen to 
be strictly attractive, and the resonant force compo- 
nents (dashed and dotted lines). With regard to the 
latter, we have separately shown the propagating part 
(q G [0,w/c) in Eq. (j8|, dashed line) and the evanes- 
cent part (q € [w/c, 00), dotted line). The rather aston- 
ishing result is that the evanescent part of the resonant 
force almost exactly cancels the non-resonant force com- 
ponent. Hence, in this highly non-equilibrium situation 
the largest contribution to the CP force arises from the 
propagating part of the resonant force. The total force 
(thick solid line in Fig. [T]) thus closely follows the latter. 
Only at very small molecule- wall separation z the force is 
given by its near-field part which, for a two-level isotropic 
molecule with hujA -C fcsT, reads 



F(r A ) 
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n(u A ) 
n(uj A ) 



|e(^)| 2 -l k B T e{0)-l 
hw A e(0) + l 



H^)+i| 2 

k B T] 
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The approximation in the second line of Eq. (|T5|) holds 
for good conductors. The force saturates in the high- 
temperature limit where the factor in square brackets ap- 



proaches — 1/2. In contrast, the non-resonant (Lifshitz- 
like) force alone would formally diverge. The predicted 
high-temperature saturation agrees with the previously 
found vanishing of the leading, linear contribution in 
k B T '/ '{hjj a) in the good-conductor limit (l5| . 

Let us next consider a molecule that is at thermal equi- 
librium with its environment, so that the probabilities p n 
are given by a Boltzmann distribution, 



Pn 



e -ftu„/(fe B T) 

. e^^/f^r) ' 



(16) 



Here, all resonant force components cancel and the force 
is given by single a non-resonant force contribution given 
by the first term in Eq. (jBJ) where the molecular polar- 
isability has to be replaced by its thermal counterpart 

m 



(17) 



In Fig. [21 we compare this equilibrium force on a ther- 
malised molecule (solid line) with the non-equilibrium 
ground-state force (dotted line) for the case of YbF. In 
contrast to LiH, both rotational (wfc„ = 9.05 x 10 10 rad/s, 
d = 1.31 x 10 -29 Cm) and vibrational transitions (uJkn — 
9.54 x 10 13 rad/s, d = 8.60 x 10~ 31 Cm) [H] give relevant 
contributions to the force, because at room temperature 
the frequency of the latter is very close to the peak of the 
spectrum uj^, n n(ujkn) determining the resonant force con- 
tributions in the retarded limit. The results for YbF are 
thus representative of those to be expected for CaF, BaF, 
LiRb, NaRb, LiCs, which also have considerable contri- 
butions from vibrational transitions at room tempera- 
ture. Figure [5] shows that in contrast to the ground-state 
force, which oscillates as a function of molecule-wall sep- 
aration (due to the influence of vibrational transitions), 
the force on a fully thermalised atom is monotonous and 
attractive (dominated by rotational transitions) . We em- 
phasise that the force at thermal equilibrium between the 
atom and its environment (solid line) is vastly overesti- 
mated by a Lifshitz-type macroscopic calculation (dashed 
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line) that uses the ground-state polarisability ao(w) as 
input parameter. The reduction factor in the near-field 
limit is approximately given by [l8[ as 



(18) 



■ Lifshitz | 



2n(wio) + 1 



for all r A . Its dependence on the relevant transition fre- 
quency clearly makes it species-dependent. The poten- 
tially very large reduction factors 1/870 for YbF at 
room temperature) imply that these molecules can be 
brought much closer to metallic surfaces than previously 
thought. 



III. DYNAMICAL CASIMIR-POLDER FORCE 

In order to understand the transition between the non- 
equilibrium ground-state force and the fully thermalised 
one, we need to investigate the full internal molecular 
dynamics in the presence of the Au surface. The time- 
dependent probabilities p n = p n (t) are governed by the 
rate equations 



Pn(t) = T nk p n {t) + ^ T knPk{t), 



where the transition rates are given by [2J 



(19) 



-J^-^lk {&(Unk)[n(0J n k) + 1] + 0{uJ kn )n(uJkn)} 



x d nk -ImG(r A ,r A , \uj nk \)-d k 



(20) 



The transition rates for LiH near an Au surface can 
easily be calculated using the Green tensor (jSJ). The re- 
sulting time-dependent probabilities p n (t) are displayed 
for the ground state and the first manifold of rotation- 
ally excited states in the lower panels in Fig. [3J with 
the respective transition matrix elements being given by 



,)/V6 



djn m_»|] m\ — du M , «0 
[24f | . For large molecule- wall separation the transition 
rates to the different substates of the first manifold are 
very similar and so are the resulting probabilities (lower 
right panel). When moving closer to the surface, the 
transition rates become affected by the evanescent and 
propagating parts of the reflected field. The contribu- 
tions of the latter are strongly oscillating so that the 
rates r|o,o>— >| i,±i> exhibit a pronounced minimum at 
z = lljum. This is not the case for the rate rio,Q>— >-|i,o) 
due to the 1/z-contribution from the evanescent fields 
(lower left panel in Fig. [3j . Hence, at first only the oc- 
cupation of the level |1,0) reaches equilibrium with the 
level 1 0, 0) , and full thermalisation is realised only at a 
much later time. 

The dynamics of the CP force is then governed by the 
internal molecular dynamics according to 



F(r A ,t) = Y,Pn(t)F n (r A ). 



(21) 
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FIG. 3: Transient CP force and internal dynamics of a LiH 
molecule initially prepared in its ground state. For details, 
see text. 



This time-dependent force is shown for a LiH molecule 
initially prepared in its ground state p n (t = 0) = 5 n o 
in the top panel of Fig. [3j We observe a gradual dis- 
appearance of the oscillating force components on a 
timescale of approximately 3s. The attractive near-field 
force reaches its equilibrium value only much later due to 
the above mentioned strongly reduced rate o>— >|i,±i> - 
Note that during the thermalisation the molecule is in 
an anisotropic state so that we have to use the general 
expression (UJ for the force components rather than its 
isotropic special case ([6]). 



IV. CONCLUSIONS AND OUTLOOK 

Studying the CP force on polar molecules near a pla- 
nar surface at finite temperature, we have found that 
even ground-state molecules are subject to resonant spa- 
tially oscillating force components at finite temperature. 
They are due to the thermal non-equilibrium between the 
molecule and its environment. A full dynamical treat- 
ment has shown that these transient forces disappear in 
the course of thermalisation of the molecule. The re- 
maining equilibrium force can be vastly different from 
that calculated using a Lifshitz-type force expression for 
ground-state molecules. 

In our numerical example of ground-state LiH, we have 
explicitly shown that the non-resonant force component 
and the evanescent part of the resonant force compo- 
nent cancel almost exactly, leaving a strongly reduced at- 
tractive force in the non-retarded limit which saturates 
at high temperatures. The force in the retarded limit 
is dominated by resonant contributions from rotational 
transitions. In contrast, the force on the heavier molecule 
YbF is dominated by resonant contributions from vibra- 



6 



tional transitions. Moreover, in thermal equilibrium at 
room temperature the resulting force is a factor 1/870 
smaller than would be expected from a Lifshitz-type cal- 
culation for the corresponding ground-state molecule. 

Whereas the CP force on a fully thermalised molecule 
is always attractive, the non-equilibrium force on a 
ground-state molecule near an Au surface at room tem- 
perature has been found to show an oscillating behaviour 
as a function of the molecule- wall separation za, with sta- 
ble equilibrium positions away from the surface. There- 
fore, one might be tempted to use these (transient) min- 
ima for trapping purposes. It turns out, however, that 
for LiH the first potential well (with its minimum at 
za = 300/im) has a depth of approximately 1CP 12 K 
which is immeasurably small. In order to increase the 
trap depth, one might envisage a situation in which the 
molecule is embedded in a planar cavity of size I con- 
sisting of two such Au surfaces. Then, the Fresnel re- 
flection coefficients in Eq. ((SJ) have to be replaced by 
r s ,p = r s _ p /(\ — r 2 sp e 2l P l ). For very good conductors such 



as Au, one can set |r SlP | « 1 — rj with t] <C 1. Hence, 
for f3l — nir (n € N), the modified Fresnel coefficients 
increase as |r ap | cx 2/rj. Choosing I close to a cavity 
resonance nir c/ to a, the contribution from propagating 
modes with small q can thus be boosted by several or- 
ders of magnitude. Thus, if by use of a cavity with a 
high Q-factor could increase the trap depth by a factor 
e.g. 10 6 , the energy difference would be in the microkelvin 
regime which could be sufficiently deep to trap cold polar 
molecules with thermal photons. This question will be 
addressed in more detail in a future investigation. 
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